Effects of Levy Flights Mobility Pattern on Epidemic Spreading under Limited 

Energy Constraint 
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Recently, many empirical studies uncovered that animal foraging, migration and human traveling 
obey Levy flights with an exponent around -2. Inspired by the deluge of H1N1 this year, in this 
paper, the effects of Levy flights' mobility pattern on epidemic spreading is studied from a network 
perspective. We construct a spatial weighted network which possesses Levy flight spatial property 
under a restriction of total energy. The energy restriction is represented by the limitation of total 
travel distance within a certain time period of an individual. We find that the exponent -2 is the 
epidemic threshold of SIS spreading dynamics. Moreover, at the threshold the speed of epidemics 
spreading is highest. The results are helpful for the understanding of the effect of mobility pattern 
on epidemic spreading. 

PACS numbers: 05.40.Fb,89.75.Ak,87.23.Ge 



I. INTRODUCTION 

With the process of globalization, the contacts between 
people of different areas become much more frequent than 
before, which makes the epidemic turn into an increas- 
ingly huge challenge for human beings. In the recent 
10 years, the worldwide deluges of epidemics happened 
in human societies have been more frequent (e.g., the 
SARS in 2003, the H5N1 in 2006 and H1N1 in 2009). 
The epidemic spreads through interactions of human or 
even animals, and it appears more powerful to damage as 
the interactions between people or animals get stronger, 
under the current global environment, the effective pre- 
caution actions should be made and taken according to 
available studies. Complex networks, as models of inter- 
actions in many real areas, such as society, technology 
and biology [l[ , have provided a practical perspective for 
studying the epidemic spreading processes accompany- 
ing the real interactions. The spreading processes can 
be regarded as dynamic processes on complex networks 
Hence, studying the characteristics and underly- 
ing mechanisms of epidemic spreading on complex net- 
works, which can be applied to a wide range of areas, 
ranging from computer virus infections Q , epidemiology 
such as the spreading of H1N1, SARS and HIV @, to 
other spreading phenomena on communication and social 
networks, such as rumor propagation Q, has attracted 
many scientists' attentions. 

There are two main aspects of such studies. One 
is to aim at setting the spreading mechanisms. For 
the epidemic spreading, the classical models include 
SIS (susccptive-infected-susceptive) model and SIR 



(susceptive- infected-remove) modelflHll Q. The other 
one is focused on researching the dynamic processes of 
epidemic spreading on networks with different topologi- 
cal structures [l|-[3| which have an essential effect on the 
dynamic processes of epidemic spreading. In regular, ran- 
dom and small- world networks, the studies of epidemic 
spreading all found that the dynamic processes undergo 
a phase transition: the effective spreading rate needs to 
exceed a critical threshold for a disease to become epi- 
demic [3, 0, Q. However, accompanied by the discovery 
that more and more networks have a scale-free distribu- 
tion of degree, the absence of a critical threshold has been 
revealed in the studies of epidemic spreading in scale- free 
networks [§4Tl|. 
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Epidemic spreading always follows the mobilities of hu- 
man and animal. Recently, many empirical studies and 
theoretical analysis on the pattern of animal foraging, mi- 
gration and human traveling have presented that these 
mobility patterns possess a Levy flights pr operty w ith an 
exponent a = — 2 (a = —1.59 for human) [L3 - fl4| . Levy 
flights means when human and animal travel, the step 
size d follows a power-law distribution P(d) ~ d a (— 3 < 
a < — 1) [l5[. Apparently, these mobility patterns can- 
not be depicted only by the topology of networks; hence, 
the characteristics of epidemic spreading following ani- 
mal and human mobility pattern should be observed on 
a network with specific spatial structure. This work is to 
aim at abstracting a network to describe the Levy flights 
mobility pattern and study the dynamic process of epi- 
demic spreading on spatial network. 

In this paper we study how the mobility pattern af- 
fects epidemic spreading from the network perspective, 
and especially, pay more attention to the extremely rapid 
spreading epidemics. We construct a weighted network 
with Levy flights spatial structure to describe the Levy 
nights mobility pattern. In this network, each node de- 
notes a small area and the weight on the edge denotes the 



2 



communication times or the quantities people or animal 
flow between the corresponding two small areas. And, as 
all the individuals only have limited energy, there must 
be a cost constraint on the mobility. Considering this 
limitation, we let the consumed energy by one commu- 
nication between any two small areas be proportional to 
the geographical distance between them. 



II. MOBILITY NETWORK 

We hypothesize that the distributions of human and 
animal mobility energy are both homogenous. Then all 
nodes can be think than they have the same energy, which 
can simplify the analysis considerably. Unfortunately, we 
have no the data about human mobility and just have 
deer and sheep moving data to support the hypothesis. 
From Fig. [IJ we can see that the distributions of con- 
sumed energy (the sum of distances for a sheep or deer 
in a day) for both sheep and deer are very narrow, thus 
the energy distribution can be looked as homogenous. 

The spatial Levy flights network can be constructed 
as follows. Based on a energy constraint uniform n 
nodes cycle (1-dimensional lattice), for each node, con- 
nections are added with power law distance distribution 
P(d) ~ d a randomly until the energy are exhausted. 
Each realization can generate a specific spatial weighted 
network. According to Levy flights mobility pattern 
(P(d) ~ d a ), the weight Wij on the link between node 
i and j should be proportional to dfj, and for a given 
network size, the sum of all Wijdij should be a constant 
which denotes the energy constraint. So, we can get an 
ensemble network model of these spatial weighted net- 
works generated by many times of realization as: 



Wij ~ d^ 

n 

V Wijd^ = W 
=i,j=i 

Solve the model we can get the network as: 



(1) 



2nE(d)J2Lid a 



(2) 



Where, W = wn denotes the total energy, w is a constant 

and E(d) = d= ,} denotes the expectation of one 

Levy flights distance. Obviously, the network is a full 
connected weighted network and each node is same with 
the degree k = ^ = 



III. 



EFFECTS OF LEVY FLIGHTS MOBILITY 
ON EPIDEMIC SPREADING 



So far we have constructed the one-dimensional Levy 
flights network. SIS epidemiological network model is 
the standard model for studying epidemic spreading Q. 



In SIS endemic network, each node has only two states, 
susceptive or infected. At each step, the susceptive node 
i is infected with rate v if connected to one infected node. 
Hence in the weighted network, the node i will be infected 
with the rate 1 — (1 — v)^i EI Wij , where I is the set of 
infected nodes. And at the next time, infected nodes are 
cured and become susceptive with rate 5. The effective 
spreading rate A, is defined as A = -|. In order to keep the 
model simple and without loss of generality, we always 
let S = 1 , which implies that all the infected nodes will be 
cured in the next step. For many networks, the epidemic 
threshold of A is a very significant index to measure the 
dynamic processes of epidemic spreading on the networks 
[E & l^-Ulf . Suppose A c is the epidemic threshold of a 
network which means that when A < A c the infection dies 
out exponentially and when A > A c the infection spreads 
and will always exist in the network. 

Suppose p denotes the fraction of infected node in the 
network, following the homogeneous mean- field approx- 
imation the dynamical rate equations for the SIS 
model arc 



^ = -p+[l-(l-A)^](l-p) 



(3) 



The first term in Eq. which is set I, is the re- 

covery rate of infected nodes. The second term takes 
into account the expected probability that a suscep- 
tive node to get the infection on this weighted network 
[I — (I — A) p - E(d > ]. Employing Taylor expansions and ignor- 
ing the infinitesimal value of higher order, when a < — 2 
the epidemic threshold A r = r = When a > — 2 

L L k w — 

and n — > +oo, E(d) — > +oo. From the Eq. (|3J), we have 
P~EW) ~* for an y P- In ordcr to keep 1-(I-A) p "eTst > 0, 
the A c must tend to 1 when n — > +oo. In conclusion we 
have 
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E(d) 



wZ(a-l) 



a > 
a < 



-2 
-2, 



(4) 



where Z(a) = 2d=i ^" i s Riemann Zeta function. Eq. 
Q shows that A c has a transition at a = —2. The 
simulated and analytical results are shown in Fig. [2] 
which are match well. More over, the results can be 
easily extended to high-dimensional space. The Eq. Q 
shows that the epidemic is liable to disappear (A c = I) 
if there are more long-distance connections than short 
ones (a > —2). This is an interesting and counterintu- 
itive conclusion. Intuitively, the epidemic will get more 
chances to spread and consequently it will infect a wider 
range of people under the condition that there are more 
long-distance connections. However, under the restric- 
tion on the total energy, more long-distance connections 
for a node mean a cost of a sharp decline in the num- 
ber of short-distance connections for exchange. Hence, 
accompanied by the decrease of the interactions between 
nodes, epidemic spreading becomes recession and liable 
to disappear. 
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FIG. 1: The frequency distribution of sum of walk distances 
in a day for deer and sheep. The data were obtained at the 
Macaulay Institute's Glensaugh Research Station, Aberdeen- 
shire, in north-east Scotland. By equipping each one of 20 
red deer and 14 sheep with a GPS collar receiver instrument, 
which is fixed to the neck of the animals, the data of locations 
during a total time period of 140 days were collected every 
30 or 60 minutes. Both distributions are very narrow and 
decay exponentially (subplots), which represents the energy 
distribution is homogenous. 
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FIG. 2: The simulated and analytical dependence of A c on a 
under w = 10. It shows the critical threshold A c obtained by 
the numeric simulation model and reaction equations under 
different exponent a. Green line shows the analytical critical 
threshold represented by Eq. (Q. As shown in Eq. (Q, there 
is a phase transition on a = —2. When a > —2, the epidemic 
threshold burst mounting up into 1, and when a < —2, the 
epidemic threshold is small values. Because when n — s> + inf , 
it is very hard to calculate A c around -2. We use a line to 
represent in the very small region around -2. 



For many epidemics such as HIV, if one is infected, 
he/she cannot be cured and always has a possibility to 
infect others. While with many other extremely rapid 
spreading epidemics such as H1N1 and SARS, we always 
ignore them or have no antibiotics to control them at 
the very beginning of spreading. In these conditions, SI 
model is a reasonable model to study the spreading pro- 



cess. There is only one difference between SIS and SI 
model. In SIS model, an infected node will have a prob- 
ability to become susceptive, while, in SI model, when a 
node is infected, it cannot be cured and always will in- 
fect other susceptive nodes. According to the outbreaks 
of H1N1 and SARS, we know that there are always a 
few infected individuals and even only one infected in- 
dividual in some cases. So, in the following numerical 
experiments, we preset only one node to be infected at 
the beginning and the effective spreading rate A = 0.05 
always. 

The spreading processes arc simulated on networks 
with different network size. At first, wc investigate the 
dependence of infected ratio on spreading time. As shown 
in Fig. [3J we can see that when a = —2, the epidemic 
spread fastest under each network size. Then, we also 
studied the dependence of terminal time T (at T, all 
nodes are infected) on a and we find that when w is not 
too large, T achieves the lowest value around the point 
of a = —2 (shown in Fig. Moreover for not too large 
energy w, T a . When a = —2, j3 « i (shown in Fig. 

ED. 

In the spatial weighted network, the spreading process 
is the most efficient when a = —2. Can wc get a similar 
conclusion in the cases of higher dimensional? It is still 
an open question because it seems tough to explain why 
epidemic spreading is fastest when a = —2 by using strict 
mathematical analysis, and the numeric simulations on 
high-dimensional networks are too expensive to do. How- 
ever, we speculate this conclusion can be extended to the 
cases of higher dimensional. Recent studies have pre- 
sented that, with the constraint W = wn, the average 
shortest path l eng th D also achieves the lowest value 
when a = — 2 |l6l. [l7j. So, wc investigate the relation- 
ship between terminal time T and the average shortest 
path length D. However, it is not convenient to obtain 
the average shortest path length for the spatial weighted 
network. The reason is that large Wij means node i and 
j are very close, that is to say, we have to transform this 
kind of weight if we want to obtain the shortest path 
which is of great importance to us. In order to avoid 
transformation, wc construct another un-wcightcd net- 
work which is similar to spatial weighted network. Based 
on a given uniform cycle network, we randomly add con- 
nections which hold the above constraint and by avoiding 
duplicated link to obtain an un-weighted network. By 
observing all phenomena detected on the above spatial 
weighted network, we find that the average shortest path 
length D indeed affects the spreading speed almost lin- 
early, shown as in Fig. [6] and Fig. [Jj although, the slops 
at the two sides of a = — 2 are slight different. This differ- 
ence implies that there are some other factors which also 
affect the spreading speed besides the average shortest 
path length. We believe that in high-dimensional net- 
works with Levy flights spatial structure, the epidemic 
also spreads fastest when a — — 2. In the future we will 
try our best to solve the question. 



4 





4000 6000 8000 10000 
Time 



FIG. 3: The dependence of spread ratio on time under dif- 
ferent a. A. Spread ratio under a = and a = —1, respec- 
tively. B. Spread ration under a = — 2. We can see that when 
a = —2, the spread ratio grows much faster than other cases. 
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FIG. 4: The dependence of terminal time T on a under 
n = 1000. It shows when the energy w is not too large, the 
terminal time T achieves the lowest value around the point 
of a = —2, while the optimal a will deviate —2 and become 
smaller when w is very large. 




FIG. 5: The dependence of terminal time T on network size 
n. It observes the impact of network size on terminal time 
under various exponent a and energy w. There is a log-log 
linear dependence between terminate time and the network 
size. Move over we find that, for not too large energy w, 
Toe re 8 . When a = -2,j9w §. 



1500 



apOOO- 




500 



FIG. 6: The relationship among the terminate time T, the 
average shortest path D and Levy flights exponent a. The 
lowest point of the curve is close at a — —2, at this point the 
average shortest path and the terminate time are both the 
minimum value. 




FIG. 7: The dependence of terminate time T on average short- 
est path D under n — 2000 with energy w = 10 and w = 20, 
the curve shown the T reach the minimum value when D 
touch the lowest point. Especially, when a > — 2, there is a 
linear relationship between D and T. 



IV. DISCUSSION AND CONCLUSION 

In summary, this paper presents the relationship be- 
tween epidemic spreading and the mobility pattern of 
human and animal from the network perspective. Ac- 
knowledging that the mobility pattern of human and an- 
imal follows the Levy nights pattern with an exponent 
a = —2, we find that the Levy nights mobility pattern is 
very efficient for epidemic spreading when a = — 2. This 
result presents a big challenge for epidemics control nowa- 
days. On the other hand, this result in a sense complies 
with the evolution theory. The special mobility pattern 
has evolved since animal appeared on the earth. Maybe 
it is very useful for searching foods, diffusing good genes 
or many others which are very important for species de- 
veloping. Consequently, it provides an effective path to 
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spread epidemic. So it is of no surprise to obtain such a 
conclusion. 
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